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periodic wave solutions are presented. And the influence of inhomogeneity coefficients on solitonic structures 
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1. Introduction 

It is important to investigate the integrability of nonlinear evolution equation (NLEE), which can be regarded 
as a pretest and the first step of its exact solvabiUty. There are many significant properties, such as bilinear 
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form, Lax pairs, infinite conservation laws, infinite symmetries, Hamiltonian structure, Painleve test and bi- 
linear Backlund transformation that can characterize integrability of nonlinear equations. Although there have 
been many methods proposed to deal with the NLEEs, e.g., inverse scattering transformation [I], Darboux 
transformation IJ], Backlund transformation(BT) ||3|, Hirota method iH) and so on. By using the bilinear form 
for a given NLEE, one can not only construct its multisoliton solutions, but also derive the bilinear BT, and 
some other properties HJ-lT]. Unfortunately, one of the key steps of this method is to replace the given NLEE 
by some more tractable bilinear equations for new Hirota's variables. There is no general rule to find the trans- 
formations, nor for choice or application of some essential formulas (such as exchange formulas). During the 
early 1930s, Bell proposed the classical Bell polynomials, which are specified by a generating function and ex- 
hibiting some important properties IS]. Since then the Bell polynomials have been exploited in combinatorics, 
statistics, and other fields ifTTI - lfTSl . However, in recent years Lambert and co-workers have proposed an al- 
ternative procedure based on the use of the Bell polynomials to obtain parameter families of bilinear Backlund 
transformation and lax pairs for soliton equations in a lucid and systematic way ISl- lfTOl . The Bell polynomials 
are found to play an important role in the characterization of integrability of a nonlinear equation. 

Recently, there has been growing interest in studying the variable-coefficient nonlinear evolution equa- 
tions (NLEEs), which are often considered to be more realistic than their constant-coefficient counterparts in 
modeUng a variety of complex nonlinear phenomena under different physical backgrounds lfT4ll . Since those 
variable-coefficient NLEEs are of practical importance, it is meaningful to systematically investigate com- 
pletely integrable properties such as bilinear form. Lax pairs, infinite conservation laws, infinite symmetries, 
Hamiltonian structure, Painleve test, bilinear Backlund transformation, symmetry algebra and construct various 
exact analytic solutions, including the soliton solutions and periodic solutions. For describing the propagation 
of solitonic waves in inhomogeneous media, the variable-coefficient KP-type equations have been derived from 
many physical apphcations in plasma physics, fluid dynamics and other fields lfT5l[T6l . 

In this paper, we will focus on a generalized variable-coefficient Kadomtsev-Petviashvili (vc-KP) equation 
with nonlinearity, dispersion and perturbed term 

[m, + hi(y, f)M3x + hiiy, t)uux]^ + hiy, t)u2x + /i4(>, t)u„ + hs(y, t)u2y + h^y, t)u,^ + hj(y, t)Uy = 0, (LI) 

where u is a differentiable function of x, y and f, hi{y, t) i - 1 , . . . , 7 are all analytic, sufficiently differentiable 
functions, may provide a more realistic model equation in several physical situations, e.g. in the propagation 
of (small-amplitude) surface waves in straits or large channels of (slowly) varying depth and width and nonva- 
nishing vorticity. Eq. ( II. lb can reduce to a series of integrable models or describe such physical phenomena 
as the electrostatic wave potential in plasma physics, the amplitude of the shallow-water wave and/or surface 
wave in fluid dynamics, etc lfT6l - |fT9l . Obviously, Eq. dl.lb contains quite a number of variable-coefficient KP 
models arising from various branches of physics, e.g. the KdV, cylindrical KdV, KP, cylindrical KP, generalized 
cylindrical KP and non-isospectral KP equations etc. Some currently important examples are given below: 
• The celebrated, historic Korteweg-de Vries (KdV) equation lfni20l 

u, + 6uuix + M3x = 0, (L2) 
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has been found to model many physical, mechanical and engineering phenomena, such as ion-acoustic waves, 
geophysical fluid dynamics, lattice dynamics, the jams in the congested traflic etc. 

• The Kadomtsev-Petviashvili (KP) equation 1211 

(m, + 6MM3.V + M3,v)j: + 0"oM2y =0, (1-3) 

where ctq = ±1, has been discovered to describe the evolution of long water waves, small-amplitude surface 
waves with weak nonlinearity, weak dispersion, and weak perturbation in the y direction, weakly relativistic 
soliton interactions in the magnetized plasma and some other nonlinear models. 

• The cyhndrical KdV equation Il22ll23]| 

u, + 6uu-ix + U3x + 'Y'^x — 0, (1-4) 

was first proposed by Maxon and Viecelli in 1974 when they studied propagation of radically ingoing acous- 
tic waves. And its counterpart in (2-i-l)-dimensional, the cylindrical KP equation Il24l l25l and generaUzed 
cylindrical KP equation IfTTl l26l 

^2 
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(u, + 6uU3y + M3v)v H ^M2v + —Ux - 0, (1-5) 

(m, + h2(t)uu2x + hi(t)uix)x + [/(O + y8(t)]u2x + r(t)u„. + — j^Mav + = 0, (1.6) 



with cTg - +1, have also been constructed to describe the nearly straight wave propagation which varies in a 
very small angular region IfTTl . Il24l - ll26l . 

• The KP equation with time-dependent coefficients ifTSl 

(m, -I- UUx + UXk)x + miOUx + I^A{f)Uly - 0, (1 .7) 

models the propagation of small-amplitude surface waves in straits or large channels of slowly varying depth 
and width and nonvanishing vorticity. 

• Jacobi elliptic function solutions and integrability property for the following variable-coefficient KP equation 

(m, + h\{t)uux + Ii2{t)u3x)t + h^(t)u2y + 6h4(t)ux — 0, (1-8) 

have been presented in Ref. 1271 . 

• The following equation 

(U, + hx{t)uUx + h2(t)U3x)x + hj,(t)U2x + h4(t)U2y - 0, (1.9) 

can be used to describe nonlinear waves with a weakly diffracted wave beam, internal waves propagating along 
the interface of two fluid layers, etc |fT9l . 

• Non-isospectral and variable-coeflScient KP equations read Il28ll 



(m, + uUx + ut,x)x + oUx + bu^. + cu2v + duxx + eu2x = 0, (1-10) 
u, + h\(uix + 6mm^ + 3cr^d^^Uyy) + h2(Ux — crxUy - Icrdl^^Uy) — h^{xUx + 2u + 2yUy) = 0, (1-1 1) 
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where a, b, c, d, e are functions of y, f, and h-, (i - 1, 2, 3) are functions of t. Bilinear representations, bilinear 
Backlund transformations and Lax pairs for non-isospectral KP equations dl.lOl ) and (II. lib are systematically 
investigated, respectively, in Refs. Il28ll . 

As we well known, the KdV, cylindrical KdV, KP, cylindrical KP, generalized cylindrical KP and non- 
isospectral KP equations belong to the integrable hierarchy of KP equation. In recent years, a large number 
of papers have been focusing on Painleve property, dromion-like structures and various exact solutions of 
NLEE Il29ll - B8ll . But their integrability, to the best of our knowledge, have not been studied in detail. The 
existence of infinite conservation laws can be considered as one of the many remarkable properties that deemed 
to characterize soliton equations. Under certain constraint conditions, the variable-coefficient models may be 
proved to be integrable and given explicit analytic solutions. The corresponding constraint conditions on Eq. 
jl.lb in this paper, which can be naturally found in the procedure of applying the Bell polynomials, will be 

h2 = CQlne^ dyli4 = h + d,lnhih2\ hs = 3a-hu dyin = dyh2 = h = 0, (1.12) 

where cq and a being both arbitrary parameters. 

The main purpose of this paper is extend the binary Bell polynomial approach to systematically construct 
bilinear formulism, bilinear Backlund transformations. Lax pairs and Darboux covariant Lax pairs of the gen- 
eralized vc-KP equation ( II. lb under conditions ( 11.121 ). To our knowledge, there have been no discussions about 
Eq. dl.lb under the conditions (11.12b . Based on its Lax equations, the infinite conservation laws of the equation 
will be constructed. By using the bilinear formula, the soliton solutions and Riemann theta function periodic 
wave solutions are also presented. 

The structure of the present paper is as follows. By virtue of the properties of the binary Bell polynomials, 
we systematically construct the biUnear representation, Backlund transformation. Lax pair and Darboux covari- 
ant Lax pairs of the generalized vc-KP equation ( 11.11 ) in Sees. 2-4, respectively. By means of its Lax equation, 
in Sec. 5, the infinite conservation laws of the equation also be constructed. In Sec. 6, based on the bilinear 
formula and the recently results in Ref. 15111521 . we present the soliton solutions and Riemann theta function 
periodic wave solutions of the generalized vc-KP equation dl.ll ) under the conditions ( 11.121 ) with co - 6. And 
we also discuss the influence of inhomogeneity coefficients on solitonic structures and interaction properties 
for physical interest and possible applications by some graphic analysis. Finally, a limiting procedure is pre- 
sented to analyze in detail, the relations between the periodic wave solutions and soliton solutions. And some 
introductions of multidimensional Bell polynomials and Riemann theta function wave are given in Appendix 
A, B, respectively. 

2. Bilinear representation 

In this section, we construct the bilinear representation of Eq. dl.lb by using an extra auxiliary variable instead 
of the exchange formulae. 
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Theorem 2.1. Using the following transformation 

M = 12/ii/i2'(ln/).„, (2.1) 

the generalized vc-KP equation dl.ll ) can be bilinearized into 

&{D,, D„ Dy) = [D,D, + lnD\ + hiDl + h4D,Dy + hsO] + (/le + d, In /zi/zj + h^^y - d]f ■ f = 0, (2.2) 

where dxf ■ f = dxf^ — ^ff^, dyf ■ f = dyf^ — 2//y, 6f ■ f = 5f^, and 6 = 6{y, t) is a constant of integration. 
Proof. To obtain the linearization of Eq. ( 11.11 ). a new variable q is introducing(^ is called a potential field) 

u = c(t)q2x, (2.3) 

where c-c{t) is a function to be determined. Substituting Eq. ( 12.3b into Eq. dl.lb . one can write the resulting 
equation of the form 

qix.t + hiqsx + ch2q2xqix + hqix + Inqix.y + h^q^oy + (/le + d, In c)^2a- + hiqxy = 0, (2.4) 

where we will see that such decomposition is necessary to get bilinear form of Eq. ( 11.11 ). Moreover by the 
integration of Eq. ( 12.4b about x, one obtains 

E{q) = qxi + hxiqAx + 3^2a.) + h^qtx + h^q^y + hsqiy + (/!& + dt In /ii/zj ')^.v + hjqy - 6, (2.5) 

by choosing the function c(f) = 6/21/12' and using the formula ( IA.7I ). where 6 = 6(y, t) is a constant of integra- 
tion. Based on the formula ( IA.7I ). Eq. ( 12.5b can be rewritten as the following form 

E(q) = Pxtiq) + hiP^xiq) + h3P2x(q) + h^Pxyiq) + hiP2y{q) + {he + d, In hih2^)qx + hjqy = 5. (2.6) 

Finally, according to the property (IA.9I ) and changing the variable 

q^2\nf^u^ c{t)q2x = Uhih^'ilnfl.x, (2.7) 

Eq. ( 12.6b produces the same bilinear representation ( 12.2b of the generalized vc-KP equation ( II. lb . □ 
The formula ( 12.21 ) is a new bilinear form, which can also reduce to the ones obtained in Refs. BllTll^Tl 
l24l I25I l49l l50l by choosing the appropriate coefficients /i, {i - 1, . . . , 7). 

(i) . If hi = (/ = 3, 4, 5, 6, 7), /11 = 1 and /12 = 6, Eq. ( 11.11 ) becomes the constant coefficient KdV equation. 
The corresponding bilinear form ( 12.2b reduces to 

[DxD, + D'jf-f^O, (2.8) 

which is also obtained in Refs. |l4l|7l|49]|50l, respectively. 

(ii) . In the case of /i, = (/ = 3,4,6,7), hi = 1, h2 = 6 and h^ = +1, Eq. ( II. lb reduces to a general KP 
equation. The corresponding bilinear form ( 12.2b becomes 

[D,D,+d4±d2]/./ = 0, (2.10) 
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which is also researched in Refs. |l4l|2T]|49], respectively. 

(iii). Assuming that hi = (/ = 3, 4, 7), /15 = 3cr^/fi and h(, = 1 /2f, Eq. ( II. lb becomes the cylindrical KP 
model II24II25I . The corresponding bilinear form ( 12.2b reduces to 



with ctq is an arbitrary constant, which is a new biUnear formulism for the cylindrical KP model. 

3. Bilinear Backlund transformation and associated Lax pair 

In this section, we construct the bilinear Backlund transformation and the Lax pair of the generalized vc- 
KP equation (11.11 ). Bilinear Backlund transformation is useful in constructing solutions and also serves as 
a characteristic of integrability for a given system. In the following, we derive a bilinear Backlund for the 
generalized vc-KP equation ( II. lb by using the use of binary Bell polynomials. 

Theorem 3.1. Suppose that f is a solution of the bilinear equation ( 12.21 ) under the conditions ( 11.12b , i.e., the 
coefficients hi (i = 1, 2, 5, 6, 7) satisfy h2 — co/Jie/'''^"'', h^ = 3a^hi, /17 = 0, then g satisfying 



is another solution of the equation (12.2b . where cq, a are arbitrary parameters and y — y(y, t) is an arbitrary 
Junction. So the sy stern (13.1b is called a bilinear Backlund transformation for the generalized vc-KP equation 



[D^D, + hiD\ + iallt^D] + (he + d, In /ij/ij ■ / = 0, 



(2.12) 



{Dl + aD,-A)f-g^O, 

[d, + hi [dI - 3aD,Dy + 3 AD,) + 1^0, + h^D, + y\f-g^Q, 



(3.1) 



(fTB . 



Proof. Suppose the following expressions 



q = 21ng, q = 21n/ 



(3.2) 



are solutions of Eq. (12.5b . respectively. The condition from the Eq. (12.5b can be changed into 



E(q') - E(q) ^(q - q),, + hi(q' - q)^, + 2>hi{q' + q)2x(q' - q)2x + hiq' - q)2x + 



iq - q)xy 



+ hsiq' - q)2y + {he + dMhxh2^){q' - q)x + h-i(q - q)y = 0. 



(3.3) 



In order to obtain such conditions, the following new auxiliary variables are introduced 



v^W- q)l^ = ln(//^), aj^W + q)l^ = ln(/g), 



(3.4) 



then we can change Eq. ( 13.3b into the following form 



E(q') - E{q) -E{(jj + v) - E{a) - v) - v_„ + h\{v4_^ + 6c02xV2x) + h^i^ix + ImVxy 



+ hsV2y + (/J6 + d, In /j 1/22 )vx + h-iVy 



=dx VS^tiv) + /zi^3.v(y, w)] + ^{v, (jj) = 0, 



(3.5) 



6 



where 

£^(v, to) = 3/ziWronskian[J%;c(i/, to), S^xiv)] + hzV2x + h^Vxy + hsV2y + (fte + d, In h\h2^)vx + h-iVy. 

To rewrite Sl{v, w) as '3^ -polynomials in form of x-divergence form and to change Eq. ( 13.51 1 into some 
conditions, one can introduce a new constant 

iViAv, w) + aiVyiv, w) = A, (3.6) 

where a = a{t) is an function of t and A is an arbitrary constant. By virtue of the Eq.( l3.6b . M{v,oj) can be 
changed into 

Si{v,(jS) - 3hiAv2x-a'^ [h5(^2x,y + (2/z5 - 3a^hi)VxVx,y + 'ia^hiV2xV^+hiV2x+ha,Vx,y+{hfy+d,\nhih2^)vx+h-iVy, 

(3.7) 

which is equivalent to the following form 

^(y, oS) = dx [OhiA + h^Wxiv) - 3ahi^x.y{v, to) + h^%{v)\ , (3.8) 

by taking 

hs - 2hs - 3a^h\ - 3a^h\, h(, + d,\nh\h2^ - 0, - 0, 

namely, 

h2 = CQhye^ hs = 3a^hu h = 0. (3.9) 
Then, using Eqs. ( I3.6b - (l3.8l l, we obtain the following system 

?V2xiv, to) + a%(v, w) - ^ = 0, 

dxiV.iv) + dx [hi [^ixiv, oj) - 3a%.(v, w) + 3A'3^,{v)] + hj^.iv) + h^^.iv)] = 0. (3.10) 

By virtue of property (IA.6I 1, Eq. (13.10b yields to the bilinear Backlund transformation (13.1b with y = y(f) is an 
arbitrary function. □ 
Backlund transformation (13.1b can be used to construct exact solutions for the generalized vc-KP equation 
dl.lb . Next, using the system ( 13.10b . we will derive Lax pairs of the equation (II. lb . 

Theorem 3.2. Under the conditions (11.12b and cq — 6, the generalized vc-KP equation dl.lb admits a Lax pair 

i^i + ady)i^ 5 i//2x + ai//y + (uef - i)i/r = 0, (3.11a) 
(d, + ^2)ip = lA/ + 4/!iiA3x - /!40'" V2.V + [6hxueS + 3hiA + h^ ^x 

+ (3/jiM.i-e/''"'" - 3hiad-^UyeS''^'" - h^a'^ ue^ + h^a'^ a) iff^O, (3.1 lb) 

where u is a solution of the equation dl.lb . 

Proof. Linearizing the Eq. d3.10b into a Lax pair, we introduce a Hopf-Cole transformation v - Using 
dA.Sb and dA.9b . one obtains 

'3^x(y) = iAa/iA, ^.vd^, w) = q2x + iAia/iA' w) = ^.t,. + lA.^./tA, 
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by means of which, Eq. ( 13.101 ) is then changed into the following form with A and y 

(^i + ady)ilf = (Aza + atAy + iqix - = 0, (3. 12a) 

{d, + S^i)^ = 4>t+ 4/2ii/'3.v - h4a'^if/2.x + (dhiqi.K + 3/ii/l + h^) if/x 

+ (3li\q3x — 3liiaq„ — hi,a^^q2x + hi,a^^ A^tjj — 0, (3.12b) 

which is equivalent to the Lax pair ( 13. Hal l and ( 13.1 Ibl ). respectively, by replacing qix with t/e/ □ 
Corollary 3.3. Using the conditions (I1.12l i and cq — 6, the Lax pair (13.1 lal l and (I3.11bb of the generalized 
vc-KP equation dl.lb is equivalent to the following Lax pair 

(ifi + ady)ifj = if/ix + ai//y + (ue/''"'" - A)i// = 0, (3.13a) 

id, + ^2)iA = >J^,- 4hiatPxy - [inue^''"" - Ih^A - /13) iAx + /t4iAy - {hiu.,e^'"^'" + Shiad'^Uye^ ''"'") (A = 0, 

(3.13b) 

where u is a solution of the equation (II. lb . 

The formulas (13.1b . ( 13.1 lab and ( 13.1 Ibb are new bilinear Backlund transformation and Lax pair, respec- 
tively, which can also reduce to the ones obtained in Refs. |[I],|I11,I1I21-|EQ1, Il24l-|l27l,|l29l,|l50l by choosing 
the appropriate coefficients /!, (/= 1, . . . , 7). Without loss of generality, taking cq - 6, then c(t) - e^f'^^'''. 

(i) . Assuming that a - hi - Q (i - 3,4,5,6,7), and hi = 1, h2 - 6, Eq. (II. lb becomes the general KdV 
model. The corresponding Backlund transformation ( 13.1b reduces to 

{Dl-A)f-g^O, 

[D,+Dl + 3ADx]f-g^0, (3.14) 

which is studied in Refs. ||4]|50l. The corresponding Lax pair ( 13.1 lab and (13.1 Ibb reduces to 

(^1 + ady)il/ = ^2x + {u- A)ip = 0, (3.15a) 
{d, + Jf2)iA = ^,+ 41A3.V + 3 (2m + A) ip, + 3m.,(A = 0, (3. 15b) 

where u is a solution of the equation dl.lb . The lax pair (I3.15ab and (I3.15bb is investigated by Lax, Ablowitz 
and co-workers in Refs. |[T]|20l, respectively. 

(ii) . For hi = (/ = 3,4,7), and hi = l/f^, /j, = h^ = 3cr2/f2, ^ Eq. dTTJ becomes the 
cylindrical KP equation Ii24ll25l . The corresponding formula (13.1b reduces to 

(D2 + croDy -A)f-g^O, 

[d, + 1 /t- {dI - 3croDxDy + 3iD,) + y] / • g = 0, (3.16) 

which is a new one and not obtained in Refs. Ii24ll25l . The corresponding Lax pair ( 13.1 lab and (13.1 Ibb reduces 
to 

i^i + ady)ip = (A2x + o-of V + (uyrt- A)>^ = 0, (3. 17a) 

id, + ^2)^ =4',+ 4/fViv + (6m ^ft/t^ + 3A/t^)if/, + (3m,. ^ft/t^ - 3cro5;'My Vf/f^)iA = 0, (3.17b) 
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where u is a solution of the equation ( 11.11 ). The lax pair ( I3.17al l and (I3.17bl l is a new one, which is not studied 
inRefs. 11241 l25l. 

(iii) . In the case of hi - II fi, /12 = /13 - f{t) + yg{t), - r{t), hs - 3crl/t^, hf, - l/2f, h^ - 0, Eq. 
jl.ll l becomes a generalized cylindrical KP equation Ifr7ll26l . The corresponding formula ( 13.1b reduces to 

[d, + [d] - 3(roD,D, + 3AD,) + (/ + yg)D, + rD, + y] / • g = 0, (3.18) 

which is also a new one and not obtained in Refs. Ifr7ll26ll . The corresponding Lax pair ( 13. Hal l and ( I3.11bb 
reduces to 

(^1 + ad,)iff = ilr2x + o-ofv + {u^t- A)>^ = 0, (3. 19a) 

(d, + ^2)iff = (A/ + 4/fV3.v - tTo'KOiAi. + [6m Vf/f' + 3A/t^ + (fit) + yg(t))] ^, 

+ [3m., ^ft/t^ - 3cro5;'Mv ^ft/t^ - o-„'r(r)M ^ft + o-oW{t)A] = 0, (3.19b) 

where m is a solution of the equation dl.lb . The lax pair (I3.19ab and (I3.19bb is a new one, which is not obtained 
inRefs. lfT7ll26l. 

(iv) . If hi = /2(f), h2 = fi(t), hs = g^{t), he = 6/(f), /i/ = (/ = 3,4,7), Eq. (O becomes a variable- 
coefficient KP equation |27|. The corresponding formula ( 13.1b reduces to 

{Dl + aiiDy -A)f-g^Q, 

[d, + 1 If- [d] - 3(roD,D, + 3AD,] + (/ + yg)D, + rD, + yj / • g = 0, (3.20) 

which is also a new one and not studied in Ref. fl]}. The corresponding Lax pair (13.1 lal ) and ( I3.11bl ) reduces 
to 

i^i + ady)if^ = if,2. + \g{t)\/ yl3m^y + (Me/'^<"'" - A)^ = 0, (3.21a) 

(d, + ^2)<A = 4/2(0«A3.v + (6/2(f)M J'-^W'" + 3/2(0^) iA. 

+ (3/2(f)M.vei' - 3/2(01^(01/ V3I(05;1mvJ -A = 0, (3.21b) 

where m is a solution of the equation dl.ll ). The lax pair ( I3.21ab and ( I3.21bb is a new one, which is not obtained 
inRefs. Il27l . 

(v) . Suppose hi - hi(t) (i - 1,2,3,5), hj = (j = 4,6,7), Eq. ( II. lb becomes a generalized variable 
coefficient KP equation lfT8l[T9ll29| . The corresponding formula ( 13.1b reduces to 

(Dl + aDy - ^)/ • ^ = 0, 

[d, + h 1 (pi - 3aD,Dy + 3 AD,) + h^D,, + 7] / ■ g = 0, (3 .22) 

which is also a new one and not obtained in Refs. lfT8l[T9ll29l . The corresponding Lax pair ( 13.1 lab and ( I3.11bb 
reduces to 

(ifi + a5v)iA = ^2x + ^lhil3hnjfy + (m - A)il/ = 0, (3.23a) 
[d, + ^2)^ = ^i+ 4/Jii/'iv + (6/iiM + 3hiA + /13) i/r., + (3/21M., - 3/21 V/i5/3/ii5;'M,.) i/^ = 0, (3.23b) 
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where u is a solution of the equation jl.ll l. The lax pair ( I3.23al l and ( I3.23bl l is a new one, which is not obtained 
inRefs. HI HH |29l . 

Starting from Lax pairs and Darboux transformation, the soliton-like solutions of the generalized vc-KP 
equation ( II. lb can be established. 

4. Darboux covariant Lax pair 

Theorem 4.1. Using the associated Lax pair ( I3.12al )-( l3.12bl ) and assuming that the parameter A is independent 
of variables x, y and t, the generalized vc-KP equation jl.ll l admits a kind of Darboux covariant Lax pair as 
follows 

+ ad,)c/> ^ A(f>, ^=dl + q2.„ (4.1a) 
{d, + ^2.cov)4> = 0, ^2.cov = 4/2i<9j - hio'^dl + {6h{q2x + /J3) + 3h{q^^ - 3hiaq„ - h4a'%^, (4. lb) 

whose fonn is Darboux covariant, namely, 

T(^i + ad,)(q)T-^ = (^+ ady)(^, (4.2a) 
Tid, + J^2.cov)iq)T-' = {d, + ^,cov)(q), (4.2b) 

with 'q — q + 2\n(p, under a certain gauge transformation 

T = 05,-0"' =5,-0-, cr = 5, In (/>. (4.3) 

The integrability condition of the Darboux covariant Lax pair (I4.1al) and (I4.1bb precisely gives rise to Eq. dl.lb 
in Lax representation 

[d, + ^2,COV, ^\ + ff'Sv] = + /!l(?4x + 3^2^) + /13?2a- + h^xy + /25?2y]x = 0, (4.4) 

if one chooses dyh/\ — h^ + dt\x\-h\hl^ , dyhi — h-j — 0. The equation ( 14.4b is equivalent to equation ( 12.5b , which 
implies that Lax equations ( 14. lab a«c/ ( I4.1bl ) is also a Lax pair for the generalized vc-KP equation dl.lb . 
Proof. Let be a solution of the Lax pair (I3.12ab . The following transformation ( 14.3b change the operator 
.^\{q) + ady - A into a new one as follows 

Ti^iiq) + ady - ^)r-' = ^(5) + ady - A, (4.5) 

which admitting the following form 

££\(q) ^ ^\(q^ q + hq), withA^ = 21n0. (4.6) 

Using transformation ( 14.3b . one should look for another one ^2,cov(^), which satisfies the following form 

^cov(?) = ^2,cov(? ^q + ^q)- (4.7) 
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Let be a solution of the following system 

(^1 + ady)^ = ifi = 5^ + q2x, (4.8a) 

id, + ^2,cov)(f> = 0, if2,cov = 4/ii5^. + /7i52 + bid, + b^, (4.8b) 

with bi (i - 1, 2, 3) are undetermined functions. To determine b, (/ =1,2, 3), one can show that ( 14.31 ) change 
d, + ^ifiOY into the following form 

T(d, + ^2.cow)T~' = 5, + ^cov, ^cov = 4/!i5^. +'£28, + h, (4.9) 

with {j = 1, 2, 3) and .ifi.cov are determined by 

'bj = bj(q) + Abj^bj(q + Aq), j ^ 1,2,3. (4.10) 

Using ( 14.31 ) and ( 14.91 ). one has 

Abi — 0, Ab2 — 12/licrv + /?! + crbi^x, 

Ai>3 = \2h\o-2x + l2h\o-crx + cri>i_j. + b2,x + 2crj:i>i. (4.1 1) 
By virtue of ( 14.101 ). one should just express bi i - 1, 2, 3 in the following form 

bj ^ Jffj(q,qx,qy,q2x,qxy,q2y,- ■ ■), y= 1,2,3, (4.12) 

and satisfies 

Aj^fj - J^jiq + Aq, qx + Aqx, qy + Aq^, ■ ■ ■ ) - Jffjiq, qx, qy, ■ ■ ■) - Abj, (4. 13) 

where Aq„^x.my = 25"' d"-\\\q,ni,n2- 1,2, . . ., and Abj can be solved by Eq. (14.11b . 
Direct calculation shows that 

b,^ci(y,t), (4.14) 

by using Eqs.( l4.1 1"I )-( I4.13I ). where c\{y, t) being an arbitrary function about y and t. 
Using Eq.( l4.13l ). one has 

Ab2 = aM. = M.,qM + M,qMx + '^2,q,My + ' ' ' = 12/iicr^- = 6hiAq2x- (4.15) 

It implies that we can determine ^2 up to an arbitrary constant C2{y, t), namely, 

h = Miqix) = 6/21^2. + C2(y, t), (4.16) 

where C2(y, t) being an arbitrary function about y and t. 
By means of Eq. (I4.8ab . one obtains 

qix = -ao-xy - (o-x + CJ^)x. (4. 17) 
Using Eas.( l4.14b . ( 14.16b and ( 14.17b into Eq.( l4.11b . one has 

Abi - €>hio-2x - bhiacTyy + 2c\0-x - ShiAq^x - 3h\aAqxy + C\Aq2x, (4.18) 
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which can be verified that the third condition 



(4.19) 



can be satisfied by choosing 



b-i = ^(^, qx, qy, qix, q.xy, qiy, ^3x, ■ ■ • ) = ^hq^x - ihxaq„ + ciiy, t)q2x + Ciiy, t). 



(4.20) 



where cjiy, t) is an arbitrary function of y and t. 

Taking ci(y, f) = -a"'/z4, C2(y, f) = /13, C3(y, f) = in Eqs.( l4.14l i, ( 14.161 ) and (I4.20l i. we obtain the Darboux 
covariant evolution equation (I4.1bl ) by using ( 14. Sab . ( I4.8bb . 

Through a tedious calculations of the Lie bracket [5, + -Sf2,cov, + ady}, one obtains the Eg. (14. 4b by 
choosing £?y/i4 = /i6 + 5, In/ii/ij', dyhi - h-j -Q. □ 

From above, we can investigate the higher ones by using the same method 

^o.cov® = 4/115';° +118"^' + ■ ■ ■ + ^, i = 5,6,7, ■■■ , (4.21) 
which can obtain other new ones of the Eq. dl.lb . 



5. Infinite conservation laws 

In this section, we derive the infinite conservation laws for the generalized vc-KP equation (II. lb by using the 
binary Bell polynomials. 

Theorem 5.1. Under the conditions (11.12b , the generalized vc-KP equation (II. lb admits an infinite conserva- 
tion laws 

^n,, + /n,x + %,y^O, «=1,2,.... (5.1) 

The conversed densities ^l^s are obtained as follows 

2^ 2 



4-- 4 
1 

'l 



1=1 



n = 2,3, 



(5.2) 



a«c/ the first fluxes ^l^s are obtained as follows 



fx = /ii =/i,2a - 6/!i(!'(9;' J^2,v + h^J'x - dhiJ'l, 



f„ = hx 



^.,2x-6 J]^A.^,+ l-*-2 2 



k=\ ki+k2+ki=n 



— 6hxa 



d-'-^n+l.y + Yj-^kd-'Jn-k.y 



+ /!3J^,, n = 3,4, 



(5.3) 
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and the second fluxes '^'^s are obtained as follows 



n 

f„ - 3/iia ^ -^k-^n-k + fi/iiQ-J^i+i + /i4^, + hsd'^^J^nj, n-2,3. 



(5.4) 



Proof. Changing (13.3b into the divergence form and using (13.5b . one can rewrite ^(u, ai) into a new form 

I%(v, (jj) - [OhiA + /i3)vv - 3hiav^Vy\x + [-3h\au)2x + h4Vx]y. (5.5) 
which is equivalent to the following form 

^^2x + + tt^v - = 0, 

+ dy ^hiav\ + h4Vx + h^Vy — 3hiaA^ - 0, (5.6) 

by using the fact dx(v,) - d,(vx) - v^t- 

Using the relationship ( 13.4b and the following new function 

ri^(q'x-qx)l2, (5.7) 

one obtains 

Ux ^7], a)x^ q.x + rj- (5.8) 
By using (15.8b into ( 15.6b . Eq. ( 13.5b can be changed into a Riccati-type equation 

qix + rix + rf- + ad:,\y- ^Q, (5.9) 

which is a new potential function about q, and a divergence-type equation 

77, + dx ('72a - '^tt' ~ ^arjd^^^riy + 6£~77) + /2377j + ^3h\ajf + /1477 + hsd'^^rjy — 3h\ae^^ - 0, (5.10) 

in which one can obtain Eq. (15.10b by virtue of the equation ( 15.9b and take A-e^. 
Introducing the following series 

?7 = £ + ^ -^niq, qx, qix, ■ ■ ■ )s'", (5.11) 

n=l 

into Eq. (15.9b and collecting the coefficients of s, one can get the formulas ( 15.2b for y„. 



In addition, substituting the expression (15.11b into Eq. (15.10b . one obtains 

\3 /CO \2 

^ .. 6e 

/7=1 V L«— 1 Vn-l 



2 ^,„s-" + 5, /ii ^ ^„,2xs-" - 2 2 ^„e-'' -6s J] ^ns-" 

11=1 [ n=l Vn=l / Vn=l 



+ 4£^^ 



+ In 



V/,= l 



-6/2 1 a 



Vn=l / V «=1 

/ 00 \2 



+ dy ^hia 
= 0, 



n=I 



Vn=l 



(5.12) 
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from which one can obtain the infinite conservation laws ( 15.11 1 

In Eq. ( 15. lb . the conversed densities are obtained by recursion formulas (I5.2l i. and the first fluxes ^l^s and 
the second fluxes '^'„s, respectively, are obtained by ( I5.3l l and ( 15.4b through a cumbersome calculation. □ 
From above, one concludes that the first fluxes ^„'s ( 15.3b and the second fluxes '^l^s ( 15.4b can be introduced 
from M, and the formula J'nj + ^„,x + ^n.v = 0, (n = 1,2,...) implies that infinite conserved densities of the 
generalized vc-KP equation ( II. lb can be obtained by using {J'n, n - 1,2,...,). Using Eqs. ( 15.2b . ( 15.3b and 
(15.4b . one can easily obtain J^„, ^„ and $^„. And the generahzed vc-KP equation dl.lb can be expressed in the 
form of the first equation for conservation law (15.1b . 



6. Soliton solution and Riemann theta function periodic wave solution 



Under the conditions ( 11.12b and cq - 6, we can discuss the solutions of the generalized vc-KP equation ( II. lb 
by using the bilinear form ( 12.2b . The following subsections are independent to each other, and the parameters 
are also independent. 



6.1 Soliton solution 



Theorem 6.1. Assuming 6-Q, under the conditions (11.12b and co = 6, the generalized vc-KP equation (II. lb 
admits a N-soliton solution as follows 



u = 12/ii/!2'(ln/)„, 

r N N 
/ = ^ exp ^p>77>+ ^ PiPjAij 



P=0,1 



Vj=i 



l<j<i<N 



(6.1) 



3/lipr/jj(pi -p -7,5 (Pi Vj -p jV,)- 



w/iere T/y = ^ijx+vjy-{hiHj+hiHj+h4Vj+h5^ij Vj)t+cj and exp(Ay) = 3;„^y(^^^^.)2_;,3(^,,^_^^^_)2 (!<;</< A^), 
w/ii'/e yu j, Vj are the parameters characterizing the j-th soliton, Yj\<j<i<N If^^ summation over all possible pairs 
chosen from N elements under the condition 1 < j < i < N, and 2p=o,i denotes the summation over all possible 
combinations of pi, pj = 0, 1 (/, j = 1 , 2, . . . , A^). 
Proof. Substituting ( 16.1b into the bilinear form ( 12.2b yields 

( N N N ' 



p=0,lp'=0,l V j=l 
( N 

X exp 



*^i=l l<j<i<N 



0, 



(6.2) 



in which the bilinear operator ^ is given by Eq.( l2.2b with 6-0. Let the coefficient of the factor 



exp 



(6.3) 
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on the left hand of i6.2\ be it follows that 

( N N N ' 

p=0,lp'=0,l V j=l j=l j=l 

' N \ 

xexp ^ {piPj+p\p])Aij =0, (6.4) 

^ 1 < j<i<N , 

where the coefficient '^{p,p') denotes that the summations over p and p' performed under the following condi- 
tions 

f 

1 - p'j, if 1 < j < m, 
p'l = 1, if m + \ < j < n, 
p\ = 0, if n + 1 < ;■ < A?. 



Pi = 



(6.5) 



By introducing a new variable 



^j^Pj-Pp 



one obtains the following equality 

' N ^ m 2 m n n n 

exp Yj (PiPj + p'iP'j)Aij = -(l+mimj)Aij + J^ Aij+ ^ ^ 

A<j<i<N J l<j<i<N 

On account of cj, , juj = +1 and the relations 



(6.6) 



(6.7) 



/=1 ,/=m+l l<j<i<N j=m+l 



^ (/ii(P? - /^J) + IniMi -f^j) + hA{vi - Vj) + hsijiJ^Vi - iJ.-.^Vj),Hj - Hi, vj - v,) 
^(-/zi(yU^ + - h^ifii - h^ivi + Vj) - h^ii^T^Vi + iJ.J^Vj),fii Vi + V;) 



exp(A,j) 



(6.8) 



one obtains 



1 



l<j<i<N 



F (/ii(//^ - yu]) + h^ifii + luivi - Vj) + hsiji. ' V/ - Hj^Vj),Hj - m, vj - v,) 

^[-hyi^l] +juj) - hiij^i+Hj) - ll4{Vi + Vj) - hsilUj^Vi + IJ.J^Vj),IUi + fij,Vi + v/j 



(6.9) 



Substituting Eqs.( l6.6b -( l6.9l l into Eg. ( 16. 4b yields 

( N N N \ 

^ ^s^Yj^ - Zl '^i^^^P^i + ^^^Pi + ''4V; + h^pTj^j), Y ^iPb ^i^i 



v=±l V !=\ 
N 



Y[ ^ {hiifJ^ - pj) + hiijdi - iij) + h^ivi - Vj) + hii/Ji^Vi - iJ.J^Vj),iJj - Vj - v,) mimj = 0, (6.10) 



where ^ - ^(exp(A,j)) is independent of the summation indices tnr, (i - 1,2, . . .,A'). If we can verify the 
identity ( 16.101 ) for ^ s 1, = 1,2,..., then ( 16.1b is the solution of Eq. ( 11.11 ). Using the bilinear form ( 12.21) . 
one can rewrite ( 16.10b as follows 



^iv(jUi, Vi,jU2, V2, . . . Viv) 



N 




' N > 


4 


' N \ 


Y ^i'^jihiii] + h^Hi + li4Vi + Ii5iij^vi)ijj - 


)- hi 


Y ^>^> 


+ h-i 


Y ^j^j 



( N \ \ N 



+I14 Y ^i^jl^iyj + hs Y ^j^j \ ri [ihiix^i/jimiiJ.i - mjHjf - hsipiVj - yU^v/)^] = 0. (6.11) 
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^Ar(//i, vi,yU2, V2, . . . , jJ-N , v^) is a Symmetric and homogeneous polynomial, and is also an even function of ^j, 
Vj {j - 1,2, . . . , A^). Suppose vi) = i+m, +V2), then we have the following relationship 

— ^ 1 — 

(6.12) 



For s 1 , « = 1,2, the identity ( 16.1 11 1 is easily verified. Let's assume that the identity hold for N -2, uti- 
lizing the relationship ( 16.12b . it is seen that ^A((yUi,yU2, . . ■ ,IUn) can be the factor by a symmetric homogeneous 
polynomial as follows 

^ N N ^ 

1=1 j<i 

(6.13) 

According to the degrees of Eqs.( l6.11l l and ( 16.131 ), ,!^N{fii,vi, . . . ,iun, v^) must be zero for £/ = 1, « > 2, and 
the identity is proved. Hence, the expression (16.11 ) is the A^-soliton solution of the generalized vc-KP equation 
(foi l. □ 

Based on the Theorem 6.1, one can easily obtain the following corollary. 
Corollary 6.2. For the case N — I, the one-soliton solution of the generalized vc-KP equation dl.lb can be 
written as follows: 

u = 12/11/22' [ln(l + e'')].„ , (6.14) 
where rj — + vy - {hi/j.^ + h^fi + h4V + h^n^^v^)t + c. For the case N — 2, the following expression 

u = 12/21/12' [ln(l +6"" +e'^'- +e'"+''-+'^'-)]^^ , (6.15) 

• 1 /7 3 J 1 J -1 2\ • f 4,, 3/li/i?/iT(/Ji-U2)"^/l5(A'lV2~/i2Vl)" 7 .7 

With T]i - UjX + v/v - (/11U7 + h^Ui + hAVj + h^u. Vf)f + c/ ; - 1,2, e '- - — — k-^, \^ , , ^— r, describes 

the two-soliton solution for equation dl.lb . 

Based on the soliton solutions obtained by the Hirota's method, we present some figures to describe the 
propagation situations of the solitary waves. Figures 1 and 2 show the pulse propagation of the fundamental 
sohton along the distance (jc, y)-surface with suitable choice of the parameters in Eg. ( 16. 14b . In Figures 3 and 4, 
we choose the same value of //i and 1U2 but different vi and V2- In this case, the phases of the two solitons are 
the same and two sets of parallel sohtons are obtained via Eq.( l6.15b . 




10 -10 
20 -20 



(a) 




-20 -10 10 20 



(b) 




(c) 




Fig. 1. (Color online) Propagation of the solitary wave for the generalized vc-KP equation i ll, lb via expression i6.14b 
with parameters: hi = l, /j2=-sech*(f), /jj = - I, hn = i, = 2, fi = \, v = 2 and c = -1. (a) Perspective view of the wave. 
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(b) Overhead view of the wave, (c) The corresponding contour plot. 




(«) (*) (c) 



Fig. 2. (Color online) Propagation of the solitary wave for the generalized vc-KP equation Jl.lt via expression J6. 14b 
with parameters: hi = y^, /i2=-sech"(f), = t, = y, hs = 2, n = I, v = 2 and c = -1. (a) Perspective view of the wave. 
(b) Overhead view of the wave, (c) The corresponding contour plot. 




(a) (i) (c) 



Fig. 3. (Color online) Evolution plots of the two solitary waves for the generalized vc-KP equation l ll.lt via expres- 
sion J6.15t with parameters: h[ = 1, A2=sech^(f), = 1, = -f, Ih = t, H\ = 1, vi = 3, /J2 = 2, V2 = 4 and C] = C2 = 0. 
(a) Perspective view of the wave, (b) Overhead view of the wave, (c) The corresponding contour plot. 




(a) (b) id) 

Fig. 4. (Color online) Evolution plots of the two solitary waves for the generalized vc-KP equation 1 11. II ) via expres- 
sion ( 16. 15t with parameters: hi = 1, /72=sech"(/), = 1, = -1, hj = t, = 1, vi = 2, /J2 = 2, vt = -2 and C[ = ci =Q. 
(a) Perspective view of the wave, (b) Overhead view of the wave, (c) The corresponding contour plot. 

6.2 Riemann theta function periodic wave solution 

Using a multidimensional Riemann theta function, in Refs. 15111521 we proposed two key theorems to systematically con- 
struct Riemann theta function periodic wave solutions for nonlinear equations and discrete soliton equations, respectively. 
Using the results in Ref. 1511 . we can directly obtain some periodic wave solutions for the generalized vc-KP equation dl.lb 
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(see details in Appendix: B). 

Considering the conditions ( I1.12t . we consider the following bilinear form when S is nonzero constant in Eg. ( 12. 21 

^(£»,, Dy, D,)f ■ f = [d,D, + hiD^ + hiDr + h.D.Dy + hsO] -s)f-f = 0. 

Let now consider the Riemann theta function 



(6.16) 



(6.17) 



where the integer value vector n = (ni,n2, . . . ,«jv)^ 6 Z'^, complex phase variables ^ = {^1,^2, ■ ■ ■ ,^nY 6 Z", and —ir is 
a positive definite and real-valued symmetric N x N matrix. 

Theorem 6.3. Assuming that i?(^, t) is a Riemann theta function for N = I with ^ = kx + ly + mt + s, the generalized vc-KP 
equation ( II. lb admits a one-periodic wave solution as follows 



where 



with 



u = l2lnK_^dl\nm^,T), 



^1^22-^2^12 e &2aii-^ia2i 
0) = , 6 = 



^11^22 - cint^ii 



^1 1^22 — ^12^21 



(6.18) 



(6.19) 



9 = e 



•\ an = Yj "12 = £ "21 = £ An\2n-\fho^' 



n=-oo n=-co 

^2 = ^ (l6/!i;r''(2;! - l)'^k'* - Ah^nHln - ifk^ - Ah^n^iln - ifkl - Ah^n^iln - iff) 
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(6.20) 



and the other parameters k, I, T and s are free. 

Proof. In order to obtain one-periodic wave solutions of Eq. dl.lt . we consider one-Riemann theta function &(^, T)asN = 1 



(6.21) 



where the phase variable ^ = kx + ly + wt + e and the parameter Imr > 0. According to the Theorem A in Appendix (see 
details in Ref. 1511 ). k, I, m and e satisfy the following system 

Y ^{Annik, Annil, 4nmm)e^"''"' = 0, (6.22a) 

/l=-oo 

Y ^(2m(2n - l)k, 2m(2n - 1)1, 2ni(2n - l)w)e<2"'-2"+'""" = 0. (6.22b) 

/j=-oo 

Substituting the bilinear form ^ | |6.16| | into system < l6.22al f. i |6.22bb yields 
+00 

Y [^^rfn^km - 256hin^7t*k^ + Wiin'TT^lc + WiAn'TT^kl + Idh^n-nH- + 6)6^"'"" = 0, (6.23a) 

n = -oo 

Y (4;r^(2'i - ifka:, - \6hin\2n - ifk"^ + Ahi7i^(2n - ifk- + 4h47r^(2n - ifkl + Ahsn-(2n - iff + 6) e(2"'-2«+i>'"^ = 0. 

n=-oo 

(6.23b) 

The notations are the same as the system ( I6.20t . the system ( I6.23at . l l6.23bt is simplified into a linear system for the 
frequency co and the integration constant 6, namely, 

\( \ ( u \ 

012 M 0\ 

(6.24) 



( 




/ "1 






flu 


O12 




w 




&i 


^ 021 


022 




, 6 j 




, b2 , 
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Now solving this system, we get a one-periodic wave solution of Eq. il.ll 

II = 1 2/11/12' 5^ In i?(^,t), 



which provided the vector (co, dY . It solves the system i6.24t with the theta function &(^, t) given by Eq.| |6.2U . The other 
parameters k, I, r and e are free. □ 
Theorem 6.4. Assuming that i?(^i,^2, t) is a Riemann theta function for N = 2 with = kiX + Ijy + (Hit + e, (; = 1, 2), the 
generalized vc-KP equation i ll. It admits a two-periodic wave solution as follows 



u = Uo + 12/j,/!2'a^lnj?(^i,f2,T), 
where the parameters coi, 0)2, Mo and 5 satisfy the linear system 

H{(i>i,ui2,uci,6)^ = b, 

with 



(6.25) 



(6.26) 



H = (h,jU^,, b = (bub2,h,hf, hn= 4;r2<2n-6»i,fc>(2ni -e,')3,(n), 

(ni.n2)eZ2 

ha= 4;r2<2n-0,-,fc)(2«2-e?)3,(n), /3,3 = - ^ 16/7i;r*<2n - 0,-, fe)-*3,(n), 

(ni,/i2)eZ2 (ni,n2)eZ^ 

hi4= ^i(n), bi= [i6h,K\2n-ei,kf -4h3n'{2n-ei,kf -4h4n^{2n-ei,k){2n-ei,l) 

-4hsn-{2n-ei,lf]5i{n), 

5i{n) = p^' • p^- ' pi^e"^", p2 = e "% pi = e " '^ i= 1,2,3,4. (6.27) 

and Oi = (6] , e]Y , 6>i = (0, 0)^, 02 = (1, 0)'^, 63 = (0, 1)^, 64 = (1, l)'^, ' = 1, 2, 3, 4, the other parameters k, /,, T,^ and e,- 
(i, 7 = 1,2) are free. 

Proof. To obtain two-periodic wave solutions of Eq. 1 11. It . we consider two-Riemann theta function &{^i,^2, t) as N = 2 



)?(^l,ft,T)= <^ 



,7Ti(T'n,,n)+27Ti{^,n) 



(6.28) 



where the phase variable ^ = (^1,^2)^^ £ = k/x + /,v + oj,f + £,, i = 1,2, n = (ni^niY 6 Z", and -/r is a positive 

definite and real-valued symmetric 2x2 matrix which can take the form 



rii ri2 

, T21 T22 , 
By considering a variable transformation 



Im(rii) > 0, Im(T22) > 0, Tiir22 - T12 < 0. 



(6.29) 



(6.30) 



u = uo + I2/21A2 a;lnj?(^i,^2,T), 
and integrating with respect to x, the becomes the following bilinear form 

^(Av, Dy, D,)f ■ f = [d,D, + hiD* + /iiWflOi + hiDl + h^D.D, + hsD; - 6) f ■ f = Q. (6.31) 
According to the Theorem B in Appendix (see details in Ref. l51l ). A:,, co, and s, (; =1,2) satisfy the following system 

2 ^(2;r(<2n - 0,-, fe>, 2mX2n - 0,-, I), 2m{2n - Oi, (rt» e'"l<^("--e;)."--e,>+<^"."->l = 0, (6.32) 
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where 0,- = {e],e']Y, Gi = (0,0Y, 02 = (l,Of, 03 = (0, if, 04 = (1, 1)'', i = 1,2,3,4. 
Substituting the bilinear form Jf i6.31b into system ( 16.32b yields 

^ [47r2<2n - 0,-, k}(2n - 0,-, w) - \(>hin\2n - 0,-, fc)* - l6hiUon\2n - 0,-, fe)"* + ^hi]x^{2n - Oi, 

+4hy-(2n - 0;, fe>(2n - 0,-, i) + 4/;5;r2(2n - 0,-, Z)^ + 5] e"[(^("-e,).".-9,>+<-rr!..7T,>] ^ = 1, 2, 3, 4. (6.33) 



The notations are the same as the system l l6.27t . Eqs. i6.33l can be written as a linear system about the frequency coi, C02, 
Uo and the integration constant 6, namely, 

hn hn hn hi4 
hi\ hxi hr>. 

/J31 /?32 /!33 '234 





Ml 












hi 




Uo 




hi 








. *4 , 



(6.34) 



Now solving this system, we get a two-periodic wave solution of Eq. dl.U 

U = Uo + 12/7i/l2'5^1n)?(^i,^2,T), 

which provided the vector (coi , 0)2, 1'o, Sf . It solves the system J6.34t with the theta function ,^2,1") given by Ea. J6.28t . 



The other parameters fc, , /, , T,y and £, (/, j = 1,2) are free. 



□ 



We now present some figures to describe the propagation situations of the periodic waves. Figure 5 shows the prop- 



agation of the one periodic wave via solution J6.18b . Figure 6 shows the propagation of the degenerate two-periodic wave 
via solution l l6.25b . And Figures 7 and 8 show the propagation of the asymmetric and symmetric two-periodic waves via 





-1 -0.5 0,5 1 



y 0,5- 




-1 -0.5 


o.'s 






-0,5- 





(a) 



(b) 



(c) 






id) 



(/) 



Fig. 5. (Color online) A one-periodic wave of the generalized vc-KP equation i ll. lb via expression l l6.18b with pa- 
rameters: hi = 1, hi = 1, /73 = 2, /i4 = 4, /is = 6, A: = 1, / = 2, T = ; and £ = 0. This figure shows that every one-periodic 
wave is one-dimensional, and it can be viewed as a superposition of overlapping solitary waves, placed one period apart, 
(a) Perspective view of the real part of the periodic wave Re(M). {b) Overhead view of the wave, the green lines are crests 
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and the red lines are troughs, (c) The corresponding contour plot, (d) Wave propagation pattern of the wave along the x 
axis, (e) Wave propagation pattern of wave along the y axis. (/) Wave propagation pattern of wave along the t axis. 




Fig. 6. (Color online) A degenerate two-periodic wave of the generalized vc-KP equation dl.lt via expression i6.25b 
with parameters: hi = 1, /j2 = 2, /ij = 4, h4 = 6, hs = S, k] = Ii = 1, fc = h = -1, Tu = i, zn = 0.5/, T22 = 2i 
and £1 = £2 = 0. This figure shows that degenerate two-periodic wave is almost one-dimensional, (a) Perspective view 
of the real part of the periodic wave Re([(). (b) Overhead view of the wave, the green points are crests and the red points 
are troughs, (c) The corresponding contour plot, (d) Wave propagation pattern of the wave along the x axis, (e) Wave 
propagation pattern of wave along the y axis. (/) Wave propagation pattern of wave along the t axis. 




-10-8-6-4-2 02 4 6 810 -10 -8 -6 -4 -2 02 4 6 810 -0.0001 -6e-05 -26-05 02e-O5 6e-05 0.0001 

X y t 



W (e) (/) 

Fig. 7. (Color online) An asymmetric two-periodic wave of the generalized vc-KP equation dl.lt via expression 
I6.25t with parameters: hj = -1, h2 = 2, = 4, = 6, = 8, fei = 0.1, /i = 1, ^2 = '2 = 0.3, Th = T12 = 0.5i, T22 = 2i 
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and El = £2 = 0. This figure shows that the asymmetric two-periodic wave is spatially periodic in three directions, but it 
need not to be periodic in either the x, y 01 1 directions, (a) Perspective view of the real part of the periodic wave Re(«). 
(b) Overhead view of the wave, the green points are crests and the red points are troughs, (c) The corresponding contour 
plot, (d) Wave propagation pattern of the wave along the x axis, (e) Wave propagation pattern of wave along the y axis. (/) 
Wave propagation pattern of wave along the t axis. 




Fig. 8. (Color online) An symmetric two-periodic wave of the generalized vc-KP equation dl.lt via expression l l6.25b 
with parameters: h[ = - \, hi = 2, = 4, ^14 = 6, A5 = 8, k[ = 1, /] = 2, ki = 3, h = 4, Th = i, T12 = 0.5i, T22 = 2; and 
£] = £2 = 0. This figure shows that the symmetric two-periodic wave is periodic in three directions, (a) Perspective view 
of the real part of the periodic wave Re(i(). (b) Overhead view of the wave, the green points are crests and the red points 
are troughs, (c) The corresponding contour plot, (d) Wave propagation pattern of the wave along the x axis, (e) Wave 
propagation pattern of wave along the y axis. (/) Wave propagation pattern of wave along the t axis. 



6.3 Asymptotic property of Riemann theta function periodic waves 

Based on the results of Ref. 1511 . the relation between the one- and two- periodic wave solutions J6. 18t , J6.25l l and the one- 
and two- soliton solutions l l6.14t . l l6.15t can be directly established as follows. 

Theorem 6.5. If the vector {u>,6)^ is a sohition of the system l l6.24t for the one-periodic wave sohition i6. 18t . we let 

k=f„ l=^„ £=™ (6.35) 
2m 2m 2m 

where fi, v and c are given in Eq. i6.14l l. TJien we have the following asymptotic properties 

5 -» 0, Ini^ Tj + JTT, i?(^, t) ^ 1 + e'', when p 0. (6.36) 

It implies that the one-periodic sohition l l6.18t converges to the one-soliton solution l l6.14t under a small amplitude limit, 
that is (m, jj) — » (k 1 , 0). 
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Proof. By using the system i6.20t . a,^, bj, i, j = 1,2, can be rewritten as the series about p 



a, i = 32n^k[p^ + + 9^^"* + • • • + n^p'"' + •••), = I + 2[p- + + p'^ + ■ ■ ■ + p^"~ +•••), 

aji = B.7T^k[p + 9p' + 25p'^ + ■■■ + {In - 1)^^^"'--"+' +••■), an = l[p^ p' + p'^ + ■ ■ ■ + p^"'-^"*^ + ■ 
bi = 32n^ [[lehin^k'* - h^k^ - h^kl - h^f) p^ + {256hiJT-k* - Ah^k- - Ah^kl - Ah^F) + ■ ■ • 
+ [lehm'^TT^k'* - hiH^k^ - h^n^kl - hsn^l^) p-"' +•••], 

b, = 8;r- [{Ahiix^k* - h^k^ - h^kl - hif) p + (32Ahi7x^k* - 9hik^ - 9h^kl - 9h^f) p^ + ■■■ 

+ (4/ii(2n - ifn-k* - hi{2n - l)V - h^{2n - ifkl - hs{2n - iff) p-"'-^"*^ +■■■]. 

With the aid of Proposition C in Appendix, we have 



(6.37) 





1 \ 
1 









f ^ ^ 

3271- k 2 




/ "1 







,00, 


, Ai = 


%n-k 2 







, As = 


^ 12n-k 2 _ 


, A3 = A4 = 0, 










32n^ii2 




1 \ 



5, = 




, B2 = 


, , 


, B5 = 


^ 72;r2A3 , 



, — ^3 — 64 — 0, 



where Ai = Ah^n^t - h^k^ - h^kl - hsP, A2 = l6hi7T^k'* - h^k^ - h^kl - h^f and A3 = 36hi7T^k* - h^k' - h^kl 
Substituting the system ( I6.38t into formulas l ID.Tt . one can obtain 



(6.38) 



- hd-. 





( r 

-r'Ai 




( 1 
8r'Ai 


1 


Xo = 




, X2 = 




, X4 = - 








32;t-Ai 


\ 



89r'Ai +9r'A3 

320;r-Ai 



Xi = X3 = 0, 



(6.39) 



From iD.2t . one then has 

CD = -r'Ai + sk-^A,p- - (89r'Ai + 9r'A3)^/ + o(p*), 

6 = 327T^Aip^ - 320;r^Ai^j'* + 0(p\ 
which implies by using relation J6.35l > that 

(5 — > 0, 2mw — > -(hi/j.^ + /?3yu + /14V + h^yT^v^), when p — » 0. 



In order to show that one-periodic wave l l6.18t degenerates to the one-soliton solution i6.14b under the limit p ■ 
expand the periodic function !?(^, t) in the form of 

)?(^, t) = 1 + (e-"^ + e-'^f) p + (e'*'"^ + e-'*'"^) p'^ + ■ ■ ■ . 

Using the transformation ( I6.35t . one has 

)?(^, t) = 1 + / + (e"? + e^) p^ + (e"-f + e^?) ja* + » 1 + when p -> 0, 

f = 2ni^ - nr = fix + vy + 2niu>t + c. 



(6.40) 

(6.41) 
0, we first 

(6.42) 



(6.43) 



Combining Eqs. ( 16.411 1 and ( 16.43b . one deduces that 

i; —> fix + vy - (hijA^ + hi/A + h^v + /j5/j"'y^)r + c, when p —> 0, 
2ni^ — > ?/ + TTT, when p — » 0. 



(6.44) 
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With the aid of Eas.l l643] l and {6A4[ . one can obtain 

&(^) 1 + e", when ^ ^ 0. (6.45) 



From above, we conclude that the one-periodic solution | |6. 18t just converges to the one-soliton solution l l6.14t as the 
amplitude go — » 0. □ 
Theorem 6.6. //(coi, eo2, wo, S)^ is a solution of the system l l6.26t for the two-periodic wave sohition ( I6.25t . we take 

ki = l^., h = i^., £; = — 7— , 1-12 = —, i = l,2, (6.46) 
2m 2m 2m 2m 

where yU„ y„ c„ / = 1, 2, and An are given in £q.l l6.15] >. Then we have the following asymptotic relations 

iiQ — > 0, (5 — > 0, 2ni£j — > + nTjj, i = 1,2, 

j?(fi,f2,T) ^ 1 + <?"' + e'"- + e"i+''2+'*n, when pupi ^ 0. (6.47) 



It implies that the two-periodic solution ( I6.25t converges to the two-soliton solution l l6.15t under a small amplitude limit, 
that is (u, pi, P2) ~* ("1j 0, 0). 

Proof. The proof issimilarto the one of Theorem 6.5. □ 

7. Conclusions and discussions 



In this paper, under the conditions I ll.l2t , we have systematically researched integrability features of the generalized vc- 
KP equation i ll. it . which is an important model of various nonlinear real situations in hydrodynamics, plasma physics 
and some other nonlinear science when the inhomogeneities of media and nonuniformities of boundaries are taken into 
consideration. Using the properties of the binary Bell polynomials, we systematically construct the bilinear representation, 
Backlund transformation. Lax pair and Darboux covariant Lax pair, respectively, which can be reduced to the ones of 
several integrable equations such as KdV l ll.2t . KP l ll.3t . cylindrical KdV l ll.4t . cylindrical KP and generalized cylindrical 
KP il.5\ equations etc. Based on its Lax equation, the infinite conservation laws of the equation also can be constructed. 
Using the bilinear formula and the recent results in Ref. 15 II 1521 . we have present the soliton solutions and Riemann theta 
function periodic wave solutions of the vc-KP equation dl.ll l. And we are also able to choose different parameters and 
functions to obtain some solutions, and also analyze their graphics in Figures 1-4 and 5-8, respectively. Finally, a limiting 
procedure is presented to analyze in detail, the relations between the periodic wave solutions and soliton solutions. In 
conclusion, the generalized vc-KP equation i ll. It is completely integrable under the conditions l ll.l2t in the sense that it 
admits bilinear Backlund transformation. Lax pair and infinite conservation laws. And the integrable constraint conditions 
( I1.12t on the variable coefficients can be naturally found in the procedure of applying binary Bell polynomials. The results 
presented in this paper may provide further evidence of structures and complete integrability of these equations. 
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Appendix A: Multidimensional Bell polynomials 

In the following, we simply recall some necessary notations on multidimensional binary Bell polynomials, for details refer, 
for instance, to Lembert and Gilson's work [8-10]. 

Suppose f=f{x\ , X2, . . . , x„) be a multi-variables function in C°°, the expression as follows 

is called muliti-dimensional Bell polynomials, where = 3^, • • • d'^^ (0 < /, < ni,i = 1,2,..., r). Taking « = 1, Bell 

polynomials are presented as follows 

Y,„if) = Y„(f,,...,f„)=y ■ • • f'" , n = y ks,,, 

>-.(/) = yiAf) = fi. + fh Y,Af) = h. + 3/.Av + /',■■■■ (A.2) 

To make the link between the Bell polynomials and the Hirota D-operator, the multi-dimensional binary Bell polyno- 
mials can be defined as follows (9) 

^„.„...*.v,.(!^,t^) = >'„„ „.(/)| I (A.3) 

h + ■■■ + 1,- is odd, 

U^/i.vi,...,;,-.v,-' /i + • • ■ + /,- is even, 

which inherit the easily recognizable partial structure of the Bell polynomials. 

To find the relationship of "3^ -polynomials and the Hirota bilinear equation ■ ■ • D"''^ F • G |4l, one should investigate 
the following identity Jg) 

^>,-v,.....«,-v, (v = InF/G, oj = InFG) = (FG)-'O'J; • ■ ■ 0'J;:F • G, (A.5) 
where F and G are both the functions of x and t. In case of F = G, Eq. JA.5t can be changed into 

10, rii + ■ ■ ■ + n,- is odd, 

(A.6) 
^«ivi.....«,Av(?). "1 +••■+'!, is even. 

By using l lA.6t and the following structure 

Plx{q) = q2x, Px.t{q) = qxt, ^4,(15') = ?4.v + Sgl,, F6.r(?) = + ISgZv^+v + 15^2.,.,.... (A.7) 

one can characterize F-polynomials. The binary Bell polynomials i%, ,:,,...,n,j,.(y, (jj) can be rewritten as P- and ^-polynomials 

(FG)-'£)'J; ■ ■ • 0'J;:F ■ G = ,...,„,.,-,- (!^,w)l.=l„f/G,<.=lnfG 

= Z Z■■■Z^("0^'^v^....'-.vX9)>'(»,-/O.^....(",■-W^^^^ (A-S) 

ni + ---+iir=even li=0 lr=0 i=Q 
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Multidimensional Bell polynomials admits the following key property 

>'n,.v,.....«,-Av(!^)li/=ln^ = l/'n, v, /lA- (A.9) 

It implies that the Hopf-Cole transformation v = \nij/, that is, if/ = F/G is a linear transformation of a-, ,...,«,.»> i^)- By 
using l lA.St and l |A.9l l, one can then construct the Lax system of the nonlinear equations. 



Appendix B: Riemann theta function periodic wave 

Based on the results in Ref. 1511 . we consider one-periodic wave solutions of nonlinear evolution equation (NLEE). Then 
Riemann theta function reduces the following Fourier series in n 

i?(^,t) = e'""'"^^""^ (B.l) 

where the phase variable § = kx] + Ixx + ■ ■ ■ + px^ + (af + e and the parameter Im(T) > 0. 

Theorem A.(Ref. l51l ) Assuming that t) is a Riemann theta function for N = \ with % = kx[ + 1x2 + ■ ■ ■ + P-^n + (at + e 
and k, I, ■ ■ ■ , p, ai, £ satisfy the following system 



2 ^ (4nmk, 4nml, ■■■ , 4nmp, 4nmw) e^"'"" = 0, 

H=-oo 

2 ^ (2m(2n - \)k, 2m{2n - 1)1, ■■■ , 2m{2n - l)p, 2m(2« - 1)cd) e'^"' 



-2n+l)OT> 



Then the following expression 



u = uo + adlln^i^). 



(B.2a) 
(B.2b) 

(B.3) 



is the one-periodic wave solution of the NLEE. 

Let us now consider the case when N=2, the Riemann theta function takes the form of 



(B.4) 



where n = (ni,n2)^ e Z^, § = (^i, ^2) 6 C^, ^, = kjXi + IjXj + • • • + p,ja; + (jo,f + £,,;= 1,2, and -ir is a positive definite 
whose real-valued symmetric 2x2 matrix is 

Til T12 

, Im(Tii) > 0, lm(T22) > 0, T,iT22 - TJj < 0. (B.5) 

[ Tn T22 , 

Theorem B. dSlI ) Assuming that i?(^i , §2, t) is one Riemann theta function with N = 2, = fe,xi + ljX2 -I- • • • -I- p,Xjv -I- u,f -I- 
E,, i = 1,2 and k,, Ij, ■ ■ ■ , p„ m,, E; (/ =1,2) satisfy the following system 

Y ^ {2m(2n - 0„ k), 2m(2n - 0„ /),■■•, 2mX2n - 0„p), 2mX2n - 0„ co)) g'"[<-^("-e,)."-e,>+<^r..ri>] ^ 

neZ- 

wtere e, = (6,', 6^)^,61 = (0,0)^, 62 = (1,0)^, 63 = (0, l)'^, 64 = (1, 1)^, i = 1,2,3,4- Then the following expression 



u = uo + aff;,ln&(h,^2). 



(C.2) 



is the two-periodic wave solution of the NLEE. 

Finally, we present a key proposition to investigate the asymptotic property of periodic waves. We write the system 



d6.24t into power series of 



flu ai2 



fl2l fl22 ) 



■ Ao + Alp +A2ip' + ■ 



(D.l) 
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- Xo + Xip + XikT + ■ 



Substituting Eqs.l lD.lt -( ID.3t into Ea. i6.24t leads to the following recursion relations 

AflXo = fio, AoX„ + A,X„-i + ■ ■ ■ + A„Xo = B,„ n>l, ne\ 

form which we then recursively get each vector X,, i = 0, 1, ■ • • . 

Proposition C. ( 1511 ) Assuming that the matrix Ao is reversible, we can obtain 



Xq - Aq Bq, X„ - Aq 
If the matrix Aq and A i are not inverse, 



B„ - ^ A,S„_ 



n > 1, ne N. 





V 


1=1 


7 














1 
















Ai = 






.00, 








2 . 



we can obtain 

where a-*'' o-*^' denote the first and second component of a two-dimensional vector a, respectively. 



(D.2) 



(D.3) 



(D.4) 



(D.5) 



(D.6) 



(D.7) 
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